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Abstract 

We demonstrate that all the known single- and double-logarithm summa- 
tions for a parton distribution function can be unified in the Collins-Soper 
resummation technique by applying soft approximations appropriate in dif- 
ferent kinematic regions to real gluon emissions. Neglecting the gluon lon- 
gitudinal momentum, we obtain the kr (double-logarithm) resummation for 
two-scale QCD processes, and the Balitsky-Fadin-Kuraev-Lipatov (single- 
logarithm) equation for one-scale processes. Neglecting the transverse mo- 
mentum, we obtain the threshold (double-logarithm) resummation for two- 
scale processes, and the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (single- 
logarithm) equation for one-scale processes. If keeping the longitudinal and 
transverse momenta simultaneously, we derive a unified resummation for 
large Bjorken variable x, and a unified evolution equation appropriate for 
both intermediate and small x. 
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1. Introduction 



It is known that radiative corrections in perturbative QCD (PQCD) pro- 
duce large logarithms at each order of the coupling contant. Double loga- 
rithms appear in processes involving two scales, such as ln 2 (p + 6) with p + the 
large longitudinal momentum of a parton and 1/6 the small inverse impact 
parameter, where b is conjugate to the parton transverse momentum kx- In 
the kinematic end-point region with large Bjorken variable x, there exist 
ln 2 (l/iV) from the Mellin transformation of ln(l — x)/(l — x) + , for which 
the two scales are the large p + and the small infrared cutoff (1 — x)p + for 
gluon emissions from the parton. Single logarithms are generated in pro- 
cesses involving only one scale, such as lnp + and ln(l/x), for which the 
relevant scales are the large p + and the small xp + , respectively. These log- 
arithmic corrections to a parton distribution function have been summed to 
all orders by various methods, which are the kr resummation for \n 2 (p + b) 
[p]], the threshold resummation for \n 2 (l/N) 0, |3|, ||, the Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi (DGLAP) equation for lnp + ||, and the Balitsky- 
Fadin-Kuraev-Lipatov (BFKL) equation for ln(l/x) ||. 

In this paper we shall demonstrate that all the above single- and double- 
logarithm summations can be derived in the Collins-Soper (CS) resummation 
technique [|l|. The point is the soft approximation for real gluon emissions, 
with which a parton distribution function <f)(x + l + /p + , |k T + 1 T |) is associ- 
ated. The arguments of <fi indicate that the parton, emerging from a hadron, 
carries the longitudinal momentum xp + + / + and the transverse momentum 
k^ + It in order to radiate a real gluon with the momentum I. If neglecting 
the l + dependence, 

<f>{x + l + /p + , |k T + 1 T |) ^ <f>(x, |k T + 1 T |) , (1) 

the scale (1 — x)p + will not appear. Hence, Eq. ([!]) is inappropriate for the 
region with large x — > 1. Under this soft approximation, we derive the kx 
resummation for intermediate x, which involves two scales: the large xp + 
and the small kx, and the BFKL equation for small x, which involves only 
one scale xp + ~ k T . If neglecting the It dependence, 

<f>(x + l + /p + , |k T + 1 T |) w <f>(x + l + /p + , k T ) , (2) 

the transverse degrees of freedom of a parton can be integrated out, and kx 
will not be a relevant scale. Therefore, Eq. (|2|) is inappropriate for small x, 
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where the scale kx is of order xp + } and not negligible. Under this soft ap- 
proximation, we derive the threshold resummation for large x, which involves 
two scales: the large xp + and the small (1 — x)p + , and the DGLAP equation 
for intermediate x, which involves only one scale xp + ~ (1 — x)p + . 

In the regions where Eqs. ([!]) and (Q) are inappropriate, we should keep 
both the l + and lx dependences, and employ <f>(x + 1 + /p + , |k T + 1 T |) for real 
gluon emissions directly In this case the three scales xp + , (1 — x)p + and kx 
exist simultaneously. We shall derive a unified resummation (a unification of 
the kx and threshold resummations) for large x, and a unified evolution equa- 
tion (a unification of the DGLAP and BFKL equations), which is suitable 
for both intermediate and small x. In conclusion, we are able to reproduce 
all the logarithmic summations and derive their unifications simply by em- 
ploying appropriate soft approximations for real gluon emissions in the CS 
technique. The results are summarized in Table I. 

2. Master Equation 

Consider a parton distribution function 4>{x, kx,p + ) for a hadron with the 
light-like momentum p^ = p + 6 fl+ , which describes the probability that a par- 
ton carries the longitudinal momentum xp + and the transverse momentum 
kx- If the parton is a quark, is written, in the minimal subtraction scheme, 
as 

<P(x,k T ,p+) = f^-J ^e-^ + y- +i ^(p\q(y-,y T )^ + q(0)\ P ) , (3) 

where 7 + is a Dirac matrix, and \p) denotes the hadron. Averages over 
spin and color are understood. If the parton is a gluon, the operator in 
the hadronic matrix element is replaced by F+(y~, yx)F^ + (0). The above 
definition is given in the axial gauge n-A = with the gauge vector n^ 1 = 
lying on the light cone. To implement the CS technique, we allow n to vary 
arbitrarily away from the light cone (n 2 ^ 0) U, and the parton distribution 
function becomes gauge dependent. However, it will be observed that the 
kernels for various logarithmic summations turn out to be n- independent. 
This is natural, since it has been shown that parton distribution functions 
defined for different n possess the same infrared structure, and thus the 
same evolutoin behavior, though different ultraviolet structure 0. After the 
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derivation, we bring n back to the light cone, and the gauge invariance of 
the parton distribution funciton is restored. That is, the arbitrary vector 
n appears only at the intermediate stage of the derivation, and acts as an 
auxiliary tool. 

The master equation in the CS technique is a differential equation of <f) in 
P + 0) HI • Because of the scale invariance of (f) in n as indicated by the gluon 
propagator, —iN^ v '(/) /I 2 , with 

N = <T ; + n 2 - — , 4 

n ■ I (n ■ I) 1 

4> depends on p + via the ratio (p ■ n) 2 /n 2 . Hence, we have the chain rule 
relating the derivative d<p/dp + to d<p/dn a , 

d n 2 d 

P + TT^ = v a —(p , (5) 

dp + v ■ n dn a 

with v a dimensionless vector along p. The operator d/dn a applies to gluon 
propagators, leading to 

d -N^ = l —{FN av + l u N" a ) . (6) 



dn a n ■ I 

The loop momentum (l u ) contracts with the vertex the differentiated gluon 
attaches, which is then replaced by a special vertex, 



v ■ nn ■ I 

This special vertex can be read off from the combination of Eqs. (|J) and (|6|). 

Employing Ward identities 0, a diagram with the contraction of Z M can 
be expressed as the difference of the diagram, in which the particle (quark or 
gluon) propagator after the contraction is removed, and the diagram, in which 
the particle propagator before the contraction is removed. Hence, a pair 
cancellation exists between the diagrams with adjacent contractions of Z M . 
The summation of all the diagrams with different differentiated gluons then 
reduces to a single new diagram, where the most external particle propagator 
is removed. That is, the special vertex appears at the outer end of the parton 
line in this new diagram. We obtain the master equation ^ , 
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shown in Fig. 1(a), where the new diagram denoted by contains the special 
vertex represented by a square. The coefficient 2 comes from the equality of 
with the special vertex on either of the two parton lines. 

The collinear region of the loop momentum I is not important because 
of the factor l/(n ■ I) in v a with nonvanishing n 2 . Therefore, the important 
regions of / are soft and hard, in which the subdiagram containing the special 
vertex can be factorized from according to Figs. 1(b) and 1(c) at lowest 
order, respectively. The second subdiagram in Fig. 1(c), as a soft subtrac- 
tion, guarantees a hard momentum flow. The remaining part is the original 
distribution function 0. Therefore, is factorized into the convolution of the 
subdiagram containing the special vertex with 0. 

The soft contribution from Fig. 1(b) is written as 



<f> s (x, k T ,p + ) = (psv(x, k T ,p + ) + (p sr (x, k T ,p + ) , 



(9) 



with 



v fi v v 



vll 2 



d 4 ~ e l 

(2k) 

X(j)(x + l + /p + ,\k T + l T \,p + ) , 



<P(x,k T ,p + ), (10) 



(11) 



corresponding to the virtual and real gluon emissions, respectively. The color 
factor C F = 4/3 should be replaced by N c = 3 in the case with the parton 
being a gluon. The additive counterterm 5K is specified in the modified 
minimal subtraction scheme. The hard contribution from Fig. 1(c) is given 
by 

(p h (x,k T ,p + ) = G{xp + 1 fi,a s (n))(i)(x,k T ,p + ) , (12) 



with the hard function 
G = -ig 2 C F fi e 



d A ~n AT & 



x 



(xp — I) 



vl 



-SG , 



as ^-C F ln XP+U 



IT 



(13) 



where 5G is an additive counterterm. In the case with the parton being a 
gluon, the expression of G can be written down straightforwardly. The gauge 
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factor v = J (v ■ n) 2 /\n 2 \ confirms our argument that depends on p + via 
the ratio (p ■ n) 2 /n 2 . 



3. kr Resummation and BFKL Equation 

We first discuss the soft approximation in Eq. (|1]) for associated with the 
real gluon emission. Fourier transforming Eq. (|Tl|) into the impact parameter 
b space in order to decouple the It integration, we derive 



4> s (x,b,p + ) = K(l/(bfi),a s (fi))4>(x,b,p A 
with the soft function 



(14) 



K 



ig 2 C F jjf 



(2vr) 4 - e 



N vP {l) 



v f3 v u 



- + 2m5{l 2 )e 



2\A\ T -h 



-5K 



ix bfi 



(15) 



Hence, in the intermediate x region involves two scales, the large xp + 
that characterizes the hard function G in Eq. (|13|) and the small 1/6 that 
characterizes the soft function K. 

Using = (f) s + cj) h , the master equation (|j) becomes 

P + ^:<fr(x,b,p + ) = 2 [K(l/(bn),a 3 (fi)) + G(xp + / fi,a s (^))~\ (p(x,b,p + ) . 

(16) 

Since both the ultraviolet divergences in K and G come from the virtual gluon 
contribution (p sv , they cancel each other, such that K + G is renormalization- 
group (RG) invariant. The single logarithms ln(6/i) and ln(xp + / fi), contained 
in K and G, respectively, are organized by the RG equations 



d K 

fi—K 



-7k 



d r 
djj, 



(17) 



The anomalous dimension of K, \x = /idSK/d/i, is given, up to two loops, 
by [£3 



lK 



a 

7Y 



s c F + 
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with rif the number of quark flavors and Ca = 3 a color factor. The solution 
of Eq. ([H]) gives 

K(l/(bfx),a a (ij)) + G(xp + /fi,a s (fi)) = - / — j K (a s (fx)) , (19) 

Jl/b H 

where the initial conditions K(l, ot a (l/6)) and G(l,a s (xp + )) that contribute 
only to the single-logarithm summation have been droped. Solving the dif- 
ferential equation (fl6|) with the above expression inserted, we obtain the kx 
resummation ||, 

<f ) (x,b,p + ) = A k (b,xp + )<j ) (0 \x) , (20) 
with the (Sudakov) exponential 

[*p + dp fP dfi 
J 1/6 P J 1/6 /i 



A k (b,xp + ) = exp 



(21) 



In the small x region with xp + ~ fey, or xp + ~ 1/6 in the 6 space, the 
above two-scale case reduces to the one-scale case. The source of double 
logarithms, i.e., the integral containing 7^ in Eq. fll9|), is less important. 
Instead of applying the RG equation (f[7D , we simply add Eqs. (|10|)-([T2]), 
or equivalently, Figs. 1(b) and 1(c). The result can be understood in the 
way that the function G introduces an ultraviolet cutoff of order xp + ~ 
kx-i which comes from the first subdiagram of Fig. 1(c), to the virtual soft 
gluon contribution. Without Fourier transformation, cf) can be reexpressed, 
according to Eqs. ( jlPD and (|TT|), as 



r d A l v f3 v l/ 
(x,k T ,p + ) = ig 2 N c J — 4 N u p(l) — 



6(k^ Zfi) , . 

■<j>\x, k T ,p + ) 



2 



+27ri5{l 2 )<f){x, |k T + l r |,p + )J , (22) 

where the color factor has been replaced by N c , because we consider the 
gluon distribution function in the small x region. The 9 function, defining 
the ultraviolet cutoff kx, is the consequence of the inclusion of G. 

To make a variation in x via a variation in p + , we assume a fixed parton 
momentum. Under this assumption, the momentum fraction x is propor- 
tional to l/p + , and we have || 

p+ ip+^ x ' kT ' p+ ^ = fcT,p+ ) • ( 23 ) 
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Performing the integrations over / + and I in Eq. (B^) and using Eq. (|23|) , 
the master equation (0) reduces to the BFKL equation JTTI, 



a. 



dcj)(x, k T ,p + 
dln(l/x) 

with the coupling constant a 



x. 



|k T + l r |, p+) - - 4)0(x, fer, , (24) 
: N c a s /n. 



4. Threshold Resummation and DGLAP Equation 

We then consider the soft approximation in Eq. (0). In this case the 
dependence on fc T can be integrated out from both sides of Eqs. (|i~0|)-(|T2|), 
and the scale (1 — x)p + enters. We employ the Mellin transformation to bring 
(p sr from the momentum fraction x space to the moment N space, 

sr (N,p + ) = [ l dxx N - 1 $ sr {x,p+) , (25) 
Jo 

under which the l + integration decouples. Combined with the soft virtual 
contribution in Eq. (|H]), we derive 

MN,P + ) = K(p + /(N^),a s (ix)) ( p(N,p + ) , (26) 

with the soft function 



1 , r d^ e l ,^v v 



K = ^j Q d Z j—N^il )v _ l 



5(1 -z) 



P 



+2m5(l 2 )6 h- z - l —j z N ~ l 



5K 



= ^Mc F ln^, (27) 

and the counterterm 5K the same as that in Eq. ([L5|). Therefore, in the large 
x region <ft involves two scales, the large xp + ~ p + from the hard function G 
in Eq. ([13]) and the small (1 — x)p + ~ p + /N from the soft function K. 

Similarly, Eqs. ([161)- ([T9|) hold but with 1/6 replaced by p + /N . To sum 
ln(l/iV), we regard the derivative p + d<p/dp + as 

dl_ p+ d4> 

P dp+ Nd{p+/N) ' 1 ; 
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which leads to the threshold resummation, 



( j>(N,p + ) = A t (iV,p + )0( o ) 

with the exponential 

A t (N,p+) =exp 



p + dp fP + dfx 
p+/N p Jp fl 



(29) 



(30) 



In the intermediate x region the above two-scale case reduces to the one- 
scale case because of xp + ~ (1 — x)p + , and the source of double logarithms 
becomes less important. Without the Mellin transformation, the addition of 
Eqs. (|i0|)-(|l2l), with the soft approximation in Eq. (H) inserted, leads to the 
DGLAP equation M, 



P 



with the kernel 



d 
dp^ 



■0(x,p H 



f fi>(*/?,p>«,p + ) 



71 



l-z) 



(31) 



(32) 



where the variable change £ = x + l + /p + has been employed. The argument 
of a s has been chosen as the single scale xp + ~ (1 — x)p + , which is of 
order p + . Note that the kernel P differs from the splitting function P qq = 
(a s C F + z 2 ) / (I — z) + by the term (z 2 — 1)/(1 — z) + , which is finite in 
the z — > 1 limit. The reason is that the real gluon emission is evaluated under 
soft approximation as deriving P, while it is calculated exactly as deriving 



P. 



qq- 



5. Unified Logarithmic Summations 

In this section we study the case, in which both the Z + and It dependences 
of (j) in Eq. (|TTD are retained. It will be shown that a unified resummation 
for large x and a unified evolution equation for intermediate and small x are 
derived. Obviously, we should apply both the Fourier and Mellin transfor- 
mations to Eq. (|TT1), and obtain 

UN,b,p + ) = K(p + /(N f i),l/(bfi),a s ( f i))4 > (N,b,p + ) , (33) 
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with the soft function 

K = ig 2 C F /j e 

+2m5{l 2 )5 
a s (fi) 



dz 



(2tt) 4 ^ 
1+ 
P 



N vP {l) 



l-z 



I 2 

-6K , 



-cv 



7T 



In— -K 



N . 



(34) 



Kq being the modified Bessel function. It is easy to examine the large b and 
N limits of the above expression: for p + b ^> N, we have K — > and the 
soft function K approaches Eq. (O) for the resummation. For N 3> p + b, 



we have K p 
resummation. 
Equation 



— \n(vp + b/N) and K approaches Eq. 



4H implies a characteristic scale of order 



iexp 



' p + b" 



for the threshold 



(35) 



Following the similar procedures from Eqs. (|T6[)-(PTD , we derive the unified 
resummation, 

( j ) (N,b,p + ) = A u (N,b,p + )^ \ (36) 

with the exponential 



A u {N,b,p + )exp 



-2 



dp 



e X p[-K {p+b/N)]/b p J exp [-K (p+b)}/b fJL 



dp, 



(37) 

which is appropriate for arbitrary b and N. The lower bound of p corresponds 
to Eq. fl3"5"|), while the lower bound of \x is motivated by the b — > oo and 6^0 
limits, at which Eq. (j37l) approaches Eq. (pi]) and Eq. (|30"D , respectively. 

In the intermediate and small x regions, it is not necessary to resum 
the double logarithms ln 2 (l/iV). After extracting the kx resummation, the 
remaining single- logarithm summation corresponds to a unification of the 
DGLAP and BFKL equations, since both the l + and It dependences have 
been kept. We reexpress the function in the integrand of <p sr , under the 
Fourier transformation, as 

(f)(x + l + /p + , b, p + ) = 0((1 - x)p + - / + )0(x, b, p + ) 

+ [(f)(x + l + /p + ,b,p + ) -6((1 -x)p + -l + )(j)(x,b J p + )} . (38) 
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The contribution from the first term is combined with <p sv , giving the soft 
function K for the kx resummation. The RG solution of K + G is given by 



K + G = a s (xp + ) ln(l - x) + \n(p + b) 



—1k{oc s {ij)) , (39) 

1/6 jJi 



where the first term on the right-hand side comes from the extra 9 function in 
Eq. fl3"8p. The color factor has been replaced by N c , since we are considering 
the gluon distribution function. The contribution from the second term is 
written as 

dH 



■N vf} {l) -2m6(l 2 )e 



2\J\ T -h 



(27r) 4 "^ w vl 
x + l + /p + , b, p + ) -6((1- x)p + - l + )4>(x, b, p + )} 

which will generate the splitting function below. 
The master equation (0) then becomes 



(40) 



P + ^0(x,b,p + ) 



x p + d/j, 

1/6 \l 



+2a s (xp~ ] 
with the splitting function 



j K (at s (ii)) - a s (xp + )\n(p + b) 
dzP gg (z)(j)(x/z,b,p + ) , 



99 



1 



'1-z) 



+ -- 2 + z(l-z) 



{x,b,p + ) 
(41) 

(42) 



obtained from Eq. fl40D . The term —2 + z(l — z) finite as z — > and z — > 1 
has been added. The term proportional to ln(l — x) in Eq. (|3~9p has been 
absorbed into P gg to arrive at the plus distribution 1/(1 — z)+. We first 
extract the exponential A from the kx resummation, 



A(x,b,Q ,p + ) =exp 1-2 J 



, xp+ dp 



I" % K (a s (»)) - a s (p)ln^\) , (43) 

Jl/b fi x\J 



ixQo P 

where Q is an arbitrary low energy scale. It is easy to justify by substitution 
that the gluon distribution function is given by 

(f>(x,b,Q) = A(x,&,Q o ,Q)0 {O) 

rQ dfi. 



+2 / dz 



-a B (xfi)A k (x, b, /I, Q)P gg (z)(f)(x/ z, 6, /i) 



(44) 
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with ft ' the initial condition of 0. Equation ( fj4| ) is the unified evolution 
equation, which can be regarded as a modified version of the Ciafaloni- 
Catani-Fiorani-Marchesini equation |T2 . 



6. Conclusion 

In this paper we have demonstrated that all the known single- and double- 
logarithm summations, including their unifications, can be derived in the CS 
technique. The point is the treatment of the real gluon contributions. Sim- 
ply adopting soft approximations appropriate in different kinematic regions, 
i.e., neglecting the l + or It dependence in the parton distribution function 
associated with the real gluon emission, the CS technique reduces to the var- 
ious logarithmic summations. If keeping both the l + and It dependences, a 
unified resummation for large x and a unified evolution equation for inter- 
mediate and small x are obtained. Our conclusion has been summarized in 
Table I. 

The kT and threshold resummations, and the DGLAP and BFKL equa- 
tions have been widely studied and applied to many QCD processes. The 
unified resummation is appropriate for the analysis of dijet production [|13|1 , 
in which the transverse energy of one jet (the trigger jet) is measured, while 
the other jet (the probe jet) has large rapidity up to 3.0, that corresponds to 
high x values. The unified evolution equation, because of its extra Q depen- 
dence at small x, is appropriate for the explanation of the HERA data of the 
proton structure function F 2 (x, Q 2 ) jnj. These subjects will be discussed 
elsewhere. 

This work was supported by the National Science Council of Republic of 
China under Grant No. NSC-88-2112-M-006-013. 
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Table I. Logarithmic summations derived from the Collins-Soper tech- 
nique under different soft approximations at different Bjorken variables x. 





small x 


intermediate x 


large x 


neglect / + 


BFKL equation 


kT resummation 




neglect It 




DGLAP equation 


threshold resummation 


no neglect 


unified 


equation 


unified resummation 



Figure Captions 

FIG. 1. (a) The derivative p + d(f>/dp + in the axial gauge, (b) The soft struc- 
ture and (c) the ultraviolet structure of the 0(a s ) subdiagram containing the 
special vertex. 
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